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Abstract 

While the dynamics for three-dimensional axially symmetric two-electron 
quantum dots with parabolic confinement potentials is in general non- 
separable we have found an exact separability with three quantum numbers 
for specific values of the magnetic field. Furthermore, it is shown that the 
magnetic properties such as the magnetic moment and the susceptibility are 
sensitive to the presence and strength of a vertical confinement. Using a 
semiclassical approach the calculation of the eigenvalues reduces to simple 
quadratures providing a transparent and almost analytical quantization of 
the quantum dot energy levels which differ from the exact energies only by a 
few percent. 

PACS numbers: 73.21. La, 03.65.Sq, 75.75. +a, 05.45.Mt 

Current nanofabrication technology allows one to control the size and shape of quantum 
dots [0-0]. Due to the confinement of the electrons in all three spatial directions the en- 
ergy spectrum is quantized creating excellent experimental and theoretical opportunities to 
study controlled single-particle and collective dynamics at the atomic scale. For example, 
depending on the experimental setup, the spectrum of a quantum dot displays shell structure 
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or follows predictions of random matrix theory (for a review see ||). Furthermore, it 
becomes possible to trace the transition from a quantum mechanical to an almost classical 
regime. 

Few-electron quantum dots have attracted special attention @|7|], since they may pro- 
vide a natural realization of a quantum bit [|I|]. The simplest quantum dot with the essential 
features of more complex systems contains two electrons. Experimental data, including 
transport measurements || and spin oscillations in the ground state under a perpendicular 
magnetic field f|, have been explained quantum mechanically as a result of the interplay 
between the two-dimensional lateral confinement potential, electron correlations and the 
magnetic field flQUnfl . While in these experiments effects of the third spatial dimension are 
somewhat hidden, they naturally come into play with a tilted magnetic field [12| . However, 



even with a perpendicular magnetic field the vertical confinement has at least two important 
consequences which will be worked out in the following: first, it changes the magnetic mo- 
ment and susceptibility with respect to the 2D results, second, the generically non-separable 
3D dynamics becomes separable for certain values of the magnetic field. We will use a 
semiclassical description which offers a simple and accurate approach to explore the effects 
of dimensionality in quantum dots. In contrast to a circular (2D) two-electron quantum 
dot whose classical dynamics is always separable and therefore regular, the corresponding 
3D-system with axial symmetry is in general a non-integrable problem with typical features 
of mixed dynamics (regular/chaotic). 

The Hamiltonian for the 3D two-electron quantum dot reads 

h = E {^(p,-^) 2 + % H^j+y") + ^ z f\ } +v G + H spm , (i) 

where Vq = ol/\^i — r 2 | is the Coulomb energy (a = e 2 / (47iee )) and H spin = g*(si + s 2 )-B 
describes the Zeeman energy. Here m* and g* are the effective electron mass and ^-factor, 
respectively, and e is the dielectric constant. The confining potential is approximated with a 
3D axially-symmetric harmonic oscillator and hu z ^ huo are the energy scales of confinement 
in the z-direction and in the xy-pl&ne, respectively. For the typical voltage ~ 1 V applied 



to the gate, the confining potential is some eV deep which is large compared to the few 
meV of the confining frequency Hence, the electron wave function is localized close 

to the minimum of the well which always can be approximated by a parabolic potential. 
In real samples the electron-electron interaction is usually screened. However, the pure 
Coulomb interaction should suffice to understand the main features of the system. For 
the perpendicular magnetic field (B || z) we choose a gauge described by the vector A = 
[Bxr]/2 = ^B(—y, x, 0). Introducing the relative and center-of-mass coordinates r = ri— T2, 
R = |(ri + r2), the Hamiltonian, Eq. ([!]), can be separated into the center-of-mass (CM) 
H CM and relative-motion (RM) H re \ terms: H = H CM + H re \ + H spin . The solution to the 
CM-Hamiltonian is well known fH| and the effect of the Zeeman energy has been discussed 
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royiTf . In the following we will concentrate on the dynamics of H re \. 
For our analysis it is convenient to use cylindrical scaled coordinates, p = p/lo, p p = 
Pplo/h, z = z/lo, p z = pJo/h, where lo = (h/pujo) 1 / 2 is the characteristic length of the 
confinement potential with the reduced mass p = m*/2. The strength parameter a of the 
Coulomb repulsion goes over to A = 2a j (HuqIo) . Using the effective mass m* = 0.067m e , the 
dielectric constant e = 12, which are typical for GaAs, and the confining frequency hujQ=3 
meV, we obtain A ~ 3. Hereafter, for the sake of simplicity, we drop the tilde, i.e. for the 
scaled variables we use the same symbols as before scaling. 

In these variables the Hamiltonian for the relative motion takes a particular simple form 
(in units of hu>o) 

H re \ 1 



e = 



fiojQ 2 
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where m = l z /h, ujl = eB/2m*c is the Larmor frequency and 

o, p =K + ^o 2 ) 1/2 (3) 

is the effective confinement frequency in the p-coordinate which depends through ujl on the 
magnetic field. 

Due to the cylindrical symmetry, the ^-component l z = p^ of the angular momentum is 
conserved and the motion in is separated from the motion in the (p, z)-plane. Since the 



Coulomb term couples the two coordinates, the problem is in general non-integrable which 
is reflected in the Poincare sections shown in Fig.l for increasing magnetic field. The chosen 



ratio uj z /luq = 3 is of the same order of magnitude as in the experiment ]T2|. For ul = and 
small values of m the motion is mainly chaotic (see Fig. la). With the magnetic field the 
frequency of oscillations along the p-coordinate can be controlled which leads to qualitatively 
different dynamical situations (Fig.lb-d). For equal effective confinement frequencies u 2 = 
u z , the Hamiltonian Eq.@ becomes separable in spherical coordinates and the dynamics is 
integrable (Fig.lc). For two other limiting cases, the dynamics is nearly integrable, namely 
in the limit m — > oo and for uj z — > oo. The latter case represents a two-dimensional quantum 
dot, classically, we have p z , z — > in this limit. 

The semiclassical quantization of the circular 2D quantum dot is particularly simple since 
it reduces to a ID WKB quantization of the p-motion due to the separability of the problem. 
For given m and p z = z = the momentum p p determined from Eq.@ enters the action 
integral 
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I P = 7T <fp P dp = - ! \p p \ dp, (4) 

ATX J 7T Jp min 

with the turning points p m i n , p ma x as the positive roots of equation p p (p) = 0. The WKB 
quantization conditions 

I p (e) = h(n p + |), n p = 0,l,..., m = 0, ±1, ... (5) 

determine the energy levels. For non- interacting electrons (A = 0) the analytical calculation 
of the action integral leads to the (quantum mechanically exact) eigen-energies 



t(2n p +|m| + l)-— m, (6) 



which are the well known Fock-Darwin energies For A ^ 0, we calculate the action 

integral Eq.(^) numerically with a few iterations to determine the quantum eigenvalues. 

The energy spectra for non-interacting and interacting electrons are shown in Fig.2. In 
the interacting case the semiclassical result, although not exact (the error is less than 1%), 
reproduces very well the quantum mechanical results |L1,10 



Turning now to the 3D quantum dot we have seen that the dynamics is separable for 
tjj\ = u 2 p = u;^ 2 + u 2 and the Hamiltonian Eq. (g) in scaled spherical coordinates takes the 
form 
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In this case the square of the total angular momentum l 2 is an additional integral of motion. 
Therefore, the classical dynamics reduces again to a one-dimensional, radial problem. Using 
Eq.(|^) and calculating the action integral for the radial motion analogous to that in Eq.([|) 
(i.e. with r instead of p), we obtain the energy levels from the standard WKB quantization 
conditions 

I r (e)=h(n r + ±), \l\=h(l + l), 

n r} l = 0,1,..., m = 0,±1,...,±/. (8) 

Note that it is only the magnetic field which generates the spherical symmetry of the problem 
and therefore its separability leading to three good quantum numbers n r , I and m. 

In the general case of an axially symmetric 3D quantum dot we have non-integrable 
motion and a semiclassical quantization is neither straight forward nor does it give results 
which allow for a simple understanding of the dynamics. Therefore, we make use of the 
fact that in real samples the confining potential in the ^-direction is much stronger than 
in the xy-pl&ne which allows us to analyze the 3D non-integrable system with the 'removal 
of resonances' method (RRM) ||14|| . To lowest order the RRM consists of averaging the 
Hamilton function over the fastest angle of the unperturbed motion (A = 0) after rewriting 
coordinates and momenta in terms of action-angle variables ( J p , J z , 6 P , 6 Z ): 
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p 2 = — [ 2j p + \m\ - 2Jj p (j p + |m|) cos 20, , (9) 
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sin 9 Z , (10) 



uJ 7 



and p p — p, p z — z. Here, j z = J z /Ti and j p = J p /h. If u z > u p one averages over the 
angle 9 Z = u z t. As a result, the motion effectively decouples into an unperturbed motion in 



the z-coordinate governed by the potential {uo z / ujq) 2 z 2 /2 and into the relative motion in the 
p-coordinate governed by the effective potential 

V*(pJ,)AMp> + £ + ±K(-2^k), (11) 



2\u J 2p 2 irp \ lo z p 2 7 
where K is the first elliptic integral. Hence, the effective Hamiltonian reads 
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e=l7 + Kff j z . (12) 

Applying a similar procedure as in the 2D case, we calculate the action integral numerically. 
The momentum p p is determined from Eq. ([l2|) and the turning points p m i n , p max are as usual 
the (positive) roots of the equation p p (p) = 0. Finally, the WKB-quantization conditions 

I p {e) = U(n p + §), 3z = n z + \, 

n p ,n z = 0,1,2,..., m = 0,±l,±2,..., (13) 

determine the energy levels. 

Comparing the exact results for eigen-energies for the spherical case uj z /uj p = 1 we 
found good agreement even for large values of the magnetic field (Fig. 3a) although RRM 
is expected to work best for oj p /uj z < 1. Without magnetic field we have uj p /uj z = 1/3 
which means that the motion in z and p approximately decouples justifying the widely 
used 2D approximation. This is also reflected in the small difference between 2D and 3D 
results (compare Fig. 2b with Fig.3a at ujl = 0). Turning on the magnetic field increases the 
coupling of the dynamics in p and z which allows the two electrons eventually to access the 
full 3D space. As a consequence, the electrons can avoid each other more effectively and the 
Coulomb interaction has a smaller effect on the 3D spectrum than on the 2D spectrum which 
is most clearly visible for the m = energies, see Fig. 3a. We can understand this effect 
quantitatively by averaging the elliptic integral in Eq. (pTTj) over the unperturbed (A = 0) 
motion in p. It gives rise to an effective charge in the Coulomb interaction Vq ~ X e s/2p, 
where 

Ks = § [K f- — 7== , ) *, (14) 
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for n p = n z = (j p — j z = 1/2). The 3D energy quantized with this effective charge for 
the repulsion is close to the full interaction (dotted line in Fig. 3a). 

The effective charge A c fr/A as a function of uj p /uj z for different m is shown in Fig.3b. 
The maximum repulsion at u) p /u z = corresponds with u z — > oo to the 2D case. The 3D 
case without magnetic field starts for our parameters oj p /uj z = 1/3 at some value \ c s/\ < 1 
which decreases further for increasing u p /u z , i.e., increasing magnetic field. This explains 
quantitatively through the effective charge the difference of the effect of a magnetic field 
on a quantum spectrum in 2D and 3D cases. However, this difference becomes weaker for 
larger m as it is seen in Fig. 3b. 

Although the ground state as a function of the magnetic field is formed piecewise by 
levels of increasing m and alternating singlet-triplet character (see Fig. 3a) the magnetic 
properties of the ground state nevertheless reveal the dimensional difference between 2D 
and 3D. At temperature T = the dot is in the ground state and the magnetic moment 
and the magnetic susceptibility are defined by /i mag = —dE gr /dB and x — d/i mag /dB, 
respectively. Both quantities exhibit discontinuities as a function of the magnetic field due 
to the symmetry changes (with respect to m and spin). We find that these discontinuities 
shift when going from the 2D quantum dot to the 3D case as shown in Fig.4. 

By relaxing the restriction of two dimensions for a quantum dot and working in the 
physical three-dimensional space we have investigated physical examples of non-integrable 
systems close to integrability. For this situation the RRM method is naturally justified, 
since the confining frequencies in quantum dots obey the condition %oo z ^> Hujq. The WKB- 
approach provides a simple and transparent way to calculate the spectrum of the 3D two- 
electron quantum dot even for small values of the quantum numbers. We have found that 
at specific values of the magnetic field uj* l = — uJq an axially-symmetric quantum dot 
exhibits spherical symmetry and its dynamics becomes completely separable with three 
integrals of motion and three corresponding quantum numbers. We have shown that the 
confinement in the ^-direction, neglected in the 2D description of quantum dots, does have 
an influence on the spectrum and consequently also on magnetic properties of the dot. In 
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fact, by changing the confining frequency in the z-direction only slightly one can increase 
or decrease the magnetic moment and the susceptibility, i.e. one can control the magnetic 
properties of the two-electron quantum dot. 
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Figure captions: 

Fig.l. Poincare surfaces of sections z = 0, p z > of the relative motion for the axially- 
symmetric 3D two-electron quantum dot {uj z /ujq = 3, A = 3, m = 0, e = 5) in the magnetic 
field for: (a) uj l = 0, (b) uj l /uj — 2.5, (c) uj l /lu = a/8 and (d) oo l /uoq = 3.3. The section 
(c) indicates that for the corresponding value of the magnetic field the system is integrable. 

Fig. 2. The energy spectrum of a circular 2D quantum dot (in units Hcuo) as a function 
of the ratio lj l /uj for n p = and m — 0, 9 in the cases: (a) A = and (b) A = 3. 

Fig. 3. (a) The comparison between energy levels (in units Houq) of the axially-symmetric 
3D quantum dot with u z /ujq = 3 and A = 3 for n p = n z = and m — 0, ...,9 obtained 
using the RRM (full lines) and exact results for the spherical case (circles). The inset shows 
a good agreement between the RRM and the exact results. The dashed and dotted lines 
display the energy level with m = for the 2D and 3D cases with A c ff at ujq/u z = and 
1/3, respectively, (b) The dependence of the effective strength of the Coulomb interaction 
A c fj/A on the ratio u) p /uj z . 

Fig. 4. Magnetic moments /i mag (a) in the units of effective Bohr magneton ji* B = 
(m e /m*) /jl b and the magnetic susceptibility x (b) f° r the 2D (dashed lines) and 3D (full 
lines) cases as a function of the magnetic field strength (in u;L/c<Jo-units). We use the same 
parameters as in Figs. 2, 3. 
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This figure "figl.gif" is available in "gif" format from: 
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